Cooperative rectification in confined Brownian ratchets 
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We analyze the rectified motion of a Brownian particle in a confined environment. We show 
the emergence of strong cooperativity between the inherent rectification of the ratchet mechanism 
and the entropic bias of the fluctuations caused by spatial confinement. Net particle transport 
may develop even in situations where separately the ratchet and the geometric restrictions do not 
give rise to particle motion. The combined rectification effects can lead to bidirectional transport 
depending on particle size, resulting in a new route for segregation. The reported mechanism can be 
used to control transport in mesostructures and nanodevices in which particles move in a reduced 
space. 
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Brownian motors, identified in a variety of conditions 
ranging from biological [l[ to synthetic systems ex- 
tract work from thermal fluctuations in out of equilib- 
rium conditions. In particular, Brownian ratchets rectify 
thermal fluctuation due to their interaction with a peri- 
odic asymmetric potential (ratchet) in a non-equilibrium 
environment. They constitute a reference class of Brown- 
ian motors and have been widely used to understand how 
molecular 0, Q as well artificial 0-11] motors operate. 
Geometrical constraints provide an alternative means to 
rectify thermal fluctuations due to the confinement they 
impose, reducing the system capability to explore space. 
Modulations in the available explored region lead to gra- 
dients in the system effective free energy, inducing a lo- 
cal bias in its diffusion that can promote a macroscopic 
net velocity for asymmetric channel profiles [7| or due 
to applied alternating fields Q. Geometric barriers con- 
stitue a common feature at small scales; they are found 
in a variety of systems, including molecular transport 
in zeolites [9fl, ionic channels [lOj, or in microfiuidic de- 
vices 111. Il2j. where their shape explains, for example, 
the magnitude of the rectifying electric signal observed 
experimentally (l3j . 

Brownian motors usually operate in spatially restricted 
environments where thermal rectification is affected by 
the geometrical constraints. Understanding such inter- 
play is then relevant in a variety of experimental situ- 
ations ranging from micrometric systems, like microfiu- 
idic devices 11 , li| or colloids moving in optical tweezer 
arrays [3], to nanometric conditions, as realized with 
molecular motors 0] , down to the atomic scale where 
optical trapping allows to manipulate cold atoms (l5| . 

In this Letter, we will show that the interplay between 
a Brownian ratchet and the geometrical constraints it ex- 
periences strongly affects the out of equilibrium dynamics 
of small particles and promote cooperative particle trans- 
port even when neither the Brownian ratchet nor the geo- 
metrical confinement rectify on their own. We will clarify 
that such net current results from the cooperative recti- 
fication of thermal fluctuations by the Brownian ratchet 




FIG. 1: Brownian ratchet and entropic barriers. A: local bi- 
ased diffusion due to confinement, described in terms of an 
effective entropic potential, S(x). B: Two-state model for a 
Brownian ratchet: In statel particles slide along the potential 
V\ and jump with rate wi ( 2 to state2 where they diffuse until 
they jump back with rate lj2,i- The jump rates operate in 
different regions along the potential period, breaking detailed 
balance. C: A Brownian motor moving in a confined environ- 
ment will be sensitive to the free energy A\ (solid) generated 
by ratchet potential Vi (dot-dashed) and entropic potential 
generated by the channel shown in panel A (dashed). 



and the geometrical constraint and that it has a signif- 
icant signal-to-noise ratio that allows to experimentally 
test such motion on length scales comparable to a few 
ratchet periods. We will also show how the interplay be- 
tween Brownian and confinement rectification can lead 
to a significant enhancement in the rectified mean veloc- 
ity and, depending on particle size, to velocity reversal 
providing a novel mechanism for particle segregation at 
the micro- and nano-scale. 

In order to identify the generic features underlying the 
cooperative rectification provided by spatial confinement 
and Brownian ratcheting, we will analyze the dynamics 
of a single particle of radius R moving in a channel with 
variable half- width h(x), where x stands for the posi- 
tion along the channel longitudinal axis, as sketched in 
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Fig. [T] Rather than analyzing explicitly the diffusion of 
a particle in such a channel under the action of a ratchet 
potential with the appropriate boundary conditions, it 
proves insightful to take advantage of the asymmetry of 
the channel geometry and describe the particle dynamics 
in terms of its displacement along the channel longitudi- 
nal axis and include the channel boundary as an entropic 
potential S(x) = k B TS(x) = k B Thi2(h(x) - R)/R 
the particle is subject to {ks stands for the Boltzmann 
constant and T is the absolute temperature). Such an 
approximation, known as Fick- Jacobs [17H19| . is exact 
for a uniform channel while its regime of validity for gen- 
tly varying confining geometries has been explicitly elu- 
cidated [201 ] . This approach has been shown very fruitful 
to understand particle transport in a variety of confined 
systems 13, 16j | . 



To address the impact of entropic restrictions on the 
Brownian motor motion, we will analyze both a flashing 
ratchet [2l[ , a generic model for the rectified motion of 
colloidal particles, and the two-state ratchet 22|, model- 
ing the rectified motion of molecular motors along biofil- 
aments which accounts effectively for the mechanochem- 
ical coupling characteristic of molecular motors 23 1 . 

A colloidal particle subject to a periodic external po- 
tential, Vi(x) expressed in units of the thermal energy 
fcgT, behaves as a flashing ratchet when the random 
force breaks detailed balance This can be sim- 

ply achieved with a Gaussian white noise with a second 

moment amplitude g{x) = \J D{x) + Q (d x Vi(x)) 2 [2lj ]. 
where Q quantifies Brownian rectification. The particle 
density, p(x), reads 



He - — J I 

dt ~ dx~ ) 2 



dx dx 



where the dimensionless free energy Ai(x) = V\{x) — 
S(x), includes the entropic potential the particle is sub- 
ject to due to the change in the number of available states 
as the channel width varies. The channel corrugation 
induces a position-dependent diffusion coefficient, D[x) 
which depends on the channel section, h(x) 0. 

The two-state ratchet model constitutes a standard, 
simple framework to describe molecular motor motion. 
As sketched in Fig. [TJ a Brownian particle jumps be- 
tween two states, i = 1, 2, which determine under which 
potential, Vi=\^, it displaces [22]] . A choice of the jump- 
ing rates Wi2,2i that breaks detailed balance, jointly with 
an asymmetric potential Vi(x), determines the average 
molecular motor velocity vq ^= 0. The conformational 
changes of the molecular motors introduce an additional 
scale which will compete with rectification and geomet- 
rical confinement. Infinitely-processive molecular mo- 
tors remain always attached to the filament along which 
they displace and are affected by the geometrical re- 
strictions only while displacing along the filament; ac- 
cordingly, we choose channel- independent binding rate 



<^2i,p{x) = k 2 \. On the contrary, highly non-processive 
molecular motors detach frequently from the biofilament 
and diffuse away, leading to a channel- driven binding 
rate U)2i,np(%) — k^x/hix). The interaction between the 
molecular motor and the biofilament is chosen for speci- 
ficity as 

Vi(x) = V [sin(27ra) + A sin (Anx)] ,d x V 2 = (1) 

in units of ksT, where the position along the filament, 
x, is expressed in units of the ratchet period, L. A deter- 
mines the asymmetry of the ratchet potential, V\ , while 
V 2 ensures free diffusion in state i = 2. Motors jump to 
the free state only in a region of width 6 around the min- 
ima of V\ , with rate uj\ 2 = fci2 ■ Accordingly, the motor 
densities, pi,p2 along the channel follow [221 ] 

dtpi(x) + d x Ji = -u)u(x)pi(x) +CJ21 (x)p 2 (x) 
d t p 2 {x)+d x J 2 = u>i 2 {x)pi{x) - io 2 \{x)p 2 {x) 

where Ji, 2 (x) = ~D(x) (d x pi i2 (x) + pi. 2 {x)d x A li2 (x)) 
stands for the current densities in each of the two states 
in which motor displaces. Depending on the motor in- 
ternal state, two dimensionless free energies, ^1,2(2;) = 
^1,2(20 — S(x), account for the interplay between the 
biofilament interaction and the channel constraints. 

The mean particle velocity is computed from the 
numerical solution of the Fokker-Planck equations 
stated above getting for the flashing ratchets vq = 
LZ^ 1 (l — e^( L >) and for the two-state model vq = 

(Ji + J 2 )L, where Z = / Q L dx^^ J* +L ' 
p D(y)vl(y). 
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To analyze the interplay between the ratchet potential 
and the entropic constraints, we will discuss a channel 
with the same periodicity as the ratchet. In particular, 
we consider that the channel half-width obeys 

h(x) = 7 + P [sin (2ttx + A0) + A sin (47re + A</>)] (2) 

where A</> accounts for the phase difference between the 
ratchet and the entropic potentials while A quantifies the 
channel asymmetry. In turn, 7 and /3, together with 
the particle radius R, control the entropic barrier height 
because the maximum change in entropy reads A5 m = 
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h(x)„ 



where (h(x) — R) the effective half-section 



h(x) min -R 

a tracer of radius R is sensitive to. 

Although when both the channel and the ratchet are 
symmetric, A = A = 0, none of them can rectify on its 
own, the interplay of these two mechanisms leads to a net 
current whose magnitude, controlled by the phase shift 
A(j>, depends on the relative position of the ratchet along 
the channel, as shown in Fig. [5J Cooperative rectifica- 
tion emerges from the entropy-driven asymmetry of the 
hopping rates from a minimum of the effective free en- 
ergy, Ai(x), to its closest minimum. Confronting panels 
A,B with panel C in Fig. [5] shows that rectification de- 
velops only when the minimum of the free energy A\ (x) 
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FIG. 2: Rectification of a Brownian motor in a symmetric 
ratchet and symmetric corrugated channel. Panel A,B: veloc- 
ity, in units of Do/L, where Do — ksT /Qn-qR and L stands for 
the potential period, of a flashing ratchet (squares) or a pro- 
cessive (circles), non-processive (triangles) motor moving ac- 
cording to the two-state model as a function of the phase-shift 
Atj> for different values of the parameter P/y — 0.7,0.8,0.9 
(the larger the symbol size, the larger f3) and fixed 7 be- 
ing the energetic potential AVi = 4.0 and AQ = 10 for 
the flashing ratchet. Panel C: free energy profile for differ- 
ent A<f>; the dotted box is the region where U12 7^ 0, the 
vertical dashed-dotted lines mark the position of the maxima 
of Ai(x) in the absence of rectification. Panel D-E: recti- 
fied velocity as a function of the entropic barrier height upon 
variation of R (filled points) and j3 (empty points) for flashing 
ratchet (squares) or a processive (circles) , non-processive (tri- 
angles) motor moving according to the two-state model. Big- 
ger points stand for bigger A<f) being A(j> = 0.1, 0.2 for the non- 
processive two-state model and flashing ratchet, Acj> — 0.1, 0.3 
for the processive two-state model. 



is not equidistant from the free energy maxima, i.e. for 
Acf> ^ nir (n = 0, 1). While for the flashing ratchet the 
net current is always in the direction of the shortest path 
the particle has to diffuse to overcome the free energy 
barrier, the intrinsic mechanism of two-state model leads 
to more involved dynamics. The interplay between in- 
ternal motor reorganization and the geometrically vary- 
ing environment can lead to qualitatively new scenarios, 
such as the reversal in the direction of motion of proces- 
sive motors, as clearly shown in panel E of Fig. [5] This 
scenario, sensitive to the effective cross section felt by 
the motor as quantified by the parameter j3, cannot be 
obtained with flashing ratchets, or with alternative ratch- 
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FIG. 3: Rectification of a Brownian motor in an asymmet- 
ric ratchet and symmetric corrugated channel. Panel A-B: 
velocity is normalized by vo, that is the intrinsic velocity of 
a particle with the same radius R under the action of the 
ratchet in the case of a flat channel with the same volume, 
symbols and parameters values as in figure [21 C-D: velocity 
as a function of /3 (open points) or R (fileed points); points 
as in fgurtj2] A<j) = 0.2, 0.3, 0.5, 0.6 for the flashing ratchet 
and A(f> — 0.1, 0.8, 0.9 for the two-state model, the larger the 
symbols the bigger Acj>. 



eting mechanisms which do not incorporate the internal 
morphological changes of the displacing particle. 

Even though the free energy, A\{x), depends explic- 
itly on the channel section, h(x), panels D,E of Fig. [2] 
show that the dependence of rectified motion of a con- 
fined particle on the channel asymmetry, /?, (open points) 
and particle size, R (filled points), are essentially cap- 
tured when expressing the rectified velocity in terms of 
the maximum entropy difference (or channel aperture), 
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Only for very small values of AS„ 
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approaching the limit of validity of the Fick- Jacobs ap- 
proximation that prescribes a faster thermalization along 
the radial direction respect to the longitudinal convec- 
tion (25[, the details of the channel shape and particle 
size affect quantitatively the net particle velocity. There- 
fore, the relative aperture of the channel, quantified by 
AS m , captures the main dependence of the particle ve- 
locity emerging from cooperative rectification and deter- 
mines, for all the ratchet models explored, the optimal 
regime for cooperative rectified motion. 

The rectified velocities displayed in Figj2] can be of 
order Dq/L for a ratchet potential of magnitude AksT 
i.e. in experimentally achievable regimes. In fact, for 
optically-driven colloids ratchet potential amplitudes one 
order of magnitude larger than the thermal energy can be 
achieved tuning the laser beam intensity while the height 
of the energy barrier molecular motors are subject to due 
to ATP hydrolysis can be as much as ~ lOfc^T 24 1. 
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For asymmetric ratchets, A ^ 0, active transport leads 
to net rectification even for a symmetric, A = 0, cor- 
rugated channel. The particle sensitivity to the chan- 
nel shape leads now to strong mean particle velocity en- 
hancement. Moreover, panels A,B of Fig. [3] show that 
confinement allows particles to move against the under- 
lying ratchet rectification. The inset of panel C and 
panel D of Fig. [3] emphasize the non-monotonic behav- 
ior of the velocity with respect to the entropy barrier 
leading to maxima of the velocity enhancement and in- 
version. Hence, AS m allows to identify, generically, the 
optimal regime for rectified transport. Such a behav- 
ior is stronger for a flashing ratchet than for molecular 
motors, leading to velocities up to 40 times larger than 
their unconfined counterparts. Since AS m depends on 
the tracer size (filled points), panels C,D of Fig. [3] in- 
dicate a novel mechanism for particle segregation based 
on particle size; depending on the phase-shift, A<p, big- 
ger particles can move faster or slower than smaller ones 
and, by fine tuning the parameters, they can be also 
trapped or even move in opposite direction. Finally, if 
particles displace in the presence of both an asymmetric 
ratchet potential and channel corrugation (A ^ 0, A ^ 0) 
the strong enhancement in the rectified particle velocity 
is kept and both the inversion in the direction of mo- 
tion and the mechanism for particle size segregation are 
generically observed for appropriate values of the ratchet 
parameters. 

In conclusion, we have studied the cooperative recti- 
fication between geometrical constraints and Brownian 
ratchets as a new mechanism for active transport in con- 
fined environments and have shown that their interplay 
profoundly affects the net motion of small particles. We 
have clarified the physical origin of such a rectification, 
which may take place even when separately neither en- 
tropic nor ratcheting can lead to net particle motion, elu- 
cidating the role played by the biased diffusion generated 
by the geometrical constraint. The rectified velocity can 
be detectable in experimentally feasible situations and 
can be strongly enhanced by increasing the amplitude of 
the ratchet barrier in the two-state model or the multi- 
plicative noise parameter Q in the flashing ratchet. In 
the presence of ratchet rectification, the entropic con- 
straints modulates the velocity in a particle-size depen- 
dent manner that leads to regimes of maximum veloc- 
ity enhancement and velocity reversal, providing a new 
mechanism for particle segregation in confined environ- 
ments. Although cooperative rectification relies on the 
phase difference between the ratchet potential and the 
corrugated channel, the particle velocity vanishes only 
when in registry. According to Fig. [5J a probability dis- 
tribution of phase shifts, p(A<fi), will still lead to rectifi- 
cation whenever */ (A(ft 2 ) <C n, with a magnitude which 
will depend on (Atfi). For larger phase shift distributions, 
-J (A(f> 2 ) ~ 7r, cooperative rectification will survive only 
for asymmetric ratchets, with A ^ 0, and will be gener- 



ically observed when both the entropic and ratchet po- 
tentials are asymmetric, A ^ 0, A ^ 0. The cooperative 
mechanism described is robust, and can be exploited to 
control active transport of particles in ionic channels [1(3] 
or nuclear pores, confined colloids 14] or nanobead trans- 
port in microfluidic devices. 
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